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A  Study  of  Some  Electromagnetic  Problems 
Relating  to  EMP  Technology 

This  is  the  final  report  under  Air  Force  Office  of  Scientific  Research 
Grant  Go.  77-3465  bearing  the  above  title  covering  the  period  from  30  Sep¬ 
tember  1977  to  23  Februarv  1979. 

The  research  covers  two  main  topics:  the  equivalent  circuit  representa¬ 
tions  of  radiating  systems  and  certain  aspects  of  sensor  characteri sti cs . 

The  detailed  treatments  have  been  written  up  and  are  attached  as  the  appendices 
to  this  report.  The  following  is  a  brief  summary  of  the  research  results. 

For  the  equivalent  circuit  problem  the  development  and  the  testing  of  a 
numerical  technique  for  rational  functions  have  been  accomplished.  A  network 
modelling  leading  to  a  canonical  ladder  configuration  is  then  formulated. 
Finally  the  equivalent  circuit  representation  of  a  thin  biconical  antenna  is 
constructed  based  on  these  methods.  The  research  shows  that  transfer  function? 
involved  in  radiating  and  scattering  problems  can  be  represented  by  such  an 
equivalent  circuit  representation .  In  contrast  to  the  equivalent  circuit  rep¬ 
resentation  suggested  by  C.E.  Baum  [Single  Port  Equi'^lent  Circuits  fnr  An¬ 
tenna  and  Scatterers,  Interaction  Mote  295,  Air  Force  Weapons  Laboratory, 

March,  1975]  our  method  does  not  require  explicit  knowledge  of  the  poles  of 
these  transfer  functions. 

t he  second  appendix  treats  the  responses  of  a  short  dipole  and  a  small 
loop  placed  in  a  right-angle  conducting  corner.  The  work  is  intended  to 
correlate  the  open-circuit  volf.agp  of  a  censor  with  the  local  surfacp  charge 
density  or  surface  current  density  on  a  scatterer.  The  right-angle  corner 
is  a  simple  structure  for  which  an  exact  formulation  of  the  problem  can 
be  given.  The  results  show  the  exact  relationship  between  the  open  circuit 


voltage  and  the  local  surface  charge  or  surface  current  density  on  the  corner. 

The  appendix  also  contains  an  investigation  of  the  impedance  functions  of  these 
probes  taking  into  consideration  the  proximity  effect  of  the  corner. 

The  research  on  the  equivalent  circuit  representation  was  presented  at 
the  nuclear  EMP  Meeting  held  at  the  University  of  New  Mexico,  June  6-8,  1978. 

The  paper  is  entitled  "A  Network  Model  for  the  Biconical  Antenna,"  by  C.B.  Sharpe 
and  C.J.  Roussi .  Our  work  on  the  sensor  research  has  just  been  completed 
and  will  be  submitted  to  the  Air  Force  Weapons  Laboratory  for  consideration  of 
publication  as  a  note  in  the  Sensor  Series. 
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Abstract 


In  this  report  the  bir.onical  antenna  is  treated  as  a  representative 
scattering  system.  It  is  shown  that  at  its  input  terminals  the  biconical 
antenna  can  be  modeled  by  a  transmission  line  terminated  in  a  cannonical  LC 
ladder  network.  The  real  and  imaginary  parts  of  the  input  impedance  of  the 
biconical  antenna  serve  as  useful  test  functions  for  studying  the  approxima¬ 
tion  of  complex  functions  of  frequency  by  rational  functions.  An  effective 
algorithm  for  this  purpose  was  implemented  and  evaluated.  It  is  also  shown 
that  over  a  limited  domain  in  the  complex  frequency  plane  the  coles  and  zero 
of  the  system  function  can  be  recovered  via  the  rational  approximation. 
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1.  Introduction 


In  developing  equivalent  circuits  for  radiating  systems  we  have 
divided  the  problem  into  two  basic  parts:  the  development  of  a  rational 
function  approximation  technique  and  the  development  of  lumoed  network 
synthesis  procedures  appropriate  to  the  system  in  question.  Of  course, 
the  first:  part  must  serve  as  a  basis  for  the  second.  Typically  the 
rational  function  will  represent  in  analytic  form  the  transfer  admittance 
of  a  system  obtained  by  experimentally  measuring  the  amplitude  ar.d  phase 
of  the  surface  current  as  a  function  of  frequency  at  some  point  on  the 
scattering  object  with  reference  to  the  incident  electric  field  at  scr^ 
reference  plane.  In  the  general  case  all  we  can  say  about  the  transfer 
function  is  that  its  poles  must  all  lie  in  the  closed  left  half-plane. 

The  zeros  may  lie  in  either  half-plane.  Among  the  parameters  in  the 
problem  are  the  polarization  and  aspect  of  the  incident  field,  the  loca¬ 
tion  and  orientation  of  the  current  probe  on  the  object  and,  of  course, 
the  shape  of  the  scattering  body.  An  important  question  which  remains 
to  be  answered  is  the  nature  of  the  dependence  of  the  poles  and  zeros 
of  the  transfer  function  on  these  parameters. 

In  order  to  explore  the  above  question,  much  of  our  initial  effort 
has  been  devoted  to  the  development  and  testing  of  a  numerical  approxi¬ 
mation  technique  for  rational  functions.  This  technique  and  its  appli¬ 
cation  in  several  representative  approximation  problems  is  described  in 
Section  3.  As  a  preliminary  exercise,  a  network  modelling  problem  lead¬ 
ing  to  a  cannonical  ladder  configuration  was  also  investigated  and  com¬ 
pleted.  This  work  is  presented  in  the  next  section.  Both  of  the  above 
studies  were  centered  around  the  biconical  antenna.  One  reason  for  this 
approach  is  t  hat  the  input  impedance  of  the  biconical  antenna  exhibits 
many  of  the  frequency  response  characteristics  of  more  general  scattering 
structures,  and  is,  therefore,  a  useful  vehicle  for  test  purposes. 

2 .  A  Metwork  .Model  for  the  Biconical  Antenna 

The  biconical  antenna  offers  an  interesting  example  cn  which  to 
test  network  modeling  techniques  because  an  exact  analytical  expression 
for  the  input  impedance  is  available.  Tai ^  has  shown  that  the  input 
impedance  at  the  center  of  the  biconical  antenna  can  be  represented 
by  a  section  of  uniform  line  terminated  in  a  frequency  -dependent  admit¬ 
tance  Y  f  ?  )  .  This  equivalent  circuit  is  illustrated  in  Figure  1,  where 
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K  denotes  the  characteristic  impedance  of  the  line.  The  following 


was  obtained  bv  Tai  for  Y, 


Yt  = 


2L  ( 


21  •'. )  +e2 1 "  [  L  (  2  3  •  )  -L  ( 42 1 )  +  2  n  2]+e  2l“  c  [  L*  (  22 1 )  -  in  2], 
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where 
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L  ( x )  = 


1-cost 
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J  J 
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s  int 


nd  the  asterisk  denotes  the  complex  conjugate.  The  real  and  imaginary 
!  are  plotted  in  Figure  [2]  for  the  case  where  the  angle 


parts  of  Y  ( 


of  the  cone,  5-  =  .01  radians. 


The  objective  of  the  work  described  in  this  section  is  to  constru 
a  lumped  network  model  for  the  load  admittance  Y^_  of  the  equivalent  cir 
cuit  shown  in  Figure  1.  Because  represents  a  positive  real  driving 
point  admittance,  it  is  possible  in  principle  to  synthesize  a  network 
model  from  the  real  part  of  Y  alone.  From  (1),  the  asymptotic  behavio 
of  the  real  par 
1  in 


if  Y  (31)  at  low  frequencies  is  given  by 


•'.-0  Re  [Y.  (3.1)]  = 


'0 


6nK‘ 


where  2^  is  the  characteristic  impedance  of  free  space.  This  suggests 


that  we  look  for  a  network  having  an  input  admittance  Y  (j . )  such  that 


Re [ Yr  ( j u ) ] 


P( 


where  P(..“)  is  a  polynomial  of  the  form 


=  ?0  +  P2'J  +  •  ’  '  +  P4n  ° 


in 


It  will  now  be  shown  that  Y  (s)  can  always  be  realized  as  an  LC  ladder 

n  2 

terminated  in  a  resistance  if  P  (o  )  >  0,  all  o . 


Suppose  a  lossless  network  is  excited  as  shown  in  Figure  [3]  an. 


assume  that  the  scattering  coefficient  32^(s)  has  the  form 


3  2  i  (  s  )  B 


(3) 


where  B(s)  is  a  polynomial  of  degree  2 n  having  all  its  roots  in  the  LHP 


Then  if  A  ( s )  is  another  polynomial  of  degree  2n,  S^(s)  will  have  the 


form 
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C  /  \  _  -Ms) 

Ml1  B(s)  ' 

unitary  condition, 

S i ,  ( s ) S 1 ^  ( - s )  +  S7^(s)S9^(-s) 


B  (  s  )  B  (  s  )  -  A  (  s  )  A  (  -  s  )  = 


The  input  admittance  can  be  empress*: 


V  [  s  ) 


1 _ S11 (S)  =  W(s) 

1  +  S  ( s  )  U  t  s  )  ' 


=  A  (s)  +  3  ( s )  and  W(s)  =  3(s)  -  A(s).  Equation  (4)  i< 


u  ( s  :  v; ;  ■ 


.  -  s :  v;  ( s )  =  2s'  , 


(s)  and  ( 3  )  .  Using"  this  result 


]  ,;w  M  )  .  w  ( —  j  ) 
2  v  M  j ... )  U  ( -  ]  ..) 


U(j  ‘ 

n 

Thus,  the  pc  :  ynom.iai  car.  be  identified  with  the  product 


r  can  bo  shown  that  the  transmission  function  given  in  (3) 


be  realized  by  a  ladder  network  of  the  form  shown  in  Figure  [ 4 ] .  The 
number  of  independent  energy  storage  elements,  2n,  corresponds  to  the 
number  of  polos  of  S? , (s)  and  the  two  zeros  of  transmission  required 
at  s-0  are  provided  by  and  L, .  It  should  be  noted  that  this  real¬ 
ization  of  S-!1  ( s )  is  not  unique. 

The  numerical  aspect  of  the  modeling  problem  involves  the  determ.in- 

•  _  2  4  2 

at:  or.  of  tr.o  coefficients  or  P  ( )  by  curve-fitting  the  function  /  P  l .  ) 

to  the  "data"  represented  by  Re  ( Y  ]  .  This  can  most  readily  be  done 

by  f orm.inv  the  objective  function 

F  =  ‘  "  W.  : - -L - P  (  .  .  2  )  j  2  ,  (6) 

1  l-'i"  1 

P.o[Yt(-d.) 

whore  the  W.  are  arbitrary  weighting  coefficients.  In  the  results  that 
f o  1  low  v; .  1  ,  ail  i,  and  the  normalization  -.=c  was  employed.  The  ortimv 
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uency  rar.;:o  over  which  the  data  function  was  sampled  is  a  function  o: 
d  was  determined  empirically. 

Of  course  the-  advantage  of  employing  a  least-squares  objective 
tier,  as  defined  in  '  .■ )  is  that  the  minimization  of  F  with  respect  to 
p. ,  1=0/  2  ...  ,  4n,  leads  to  a  system  of  linear  equations  for  the 

own  coefficients .  The  algorithm,  is  easily  imolementei  and  the  solu- 


;f f  icients  presents  no  difficulty  as  long  as  the  corres- 
rem.ains  wel  1-conditioned .  Chebyshev  polynomial  meth.cus 
approximations  such  as  Lawson's  algorithm.  [2],  can  be 


or  r. ,  aJLtnopon  v;e  a: 


■  that  frequently  occurs  for  lar: 
in d  this  to  be  necessary .  Once 


?  [  .  )  is  known  U ( s )  can  be  recovered  from 


)  =  U ( s ) U  (-s' 


■J(s;  ha: 


van  nave  tne  norm 


=  C  .  S  1-  c  _  s 


2n-l 


2n-l‘ 


it  ter  polynomial  car.  be  determined  uniquely  from  U(s)  by  imposing 

onditicn  given  in  (5).  For  instance,  in  the  case  of  n=3,  the 

ions  for  the  c.  take  the  form, 
i 


sis  procedure  is  completed  by  expanding  Y  (s)  =  W (s ) /U ( s )  in 

n 

:  fraction  and  identifying  the  coefficients  with  the  elements 
d-,-r  network.  Thus,  for  the  network  of  Figure  4,  the  expansion 
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Figures  5  and  6  illustrate  some  typical  results  obtained  by  the 

modeling  technique  described  above.  In  each  case  9 q  =  .01  radians.  The 

element  values  for  the  approximating  network  are  calculated  assuming 

>',=  lm.  It  can  be  seen  that  the  least-squares  fit  in  the  real  part  of  the 

acm: ttance  is  quite  satisfactory  and  that  the  band  over  which  the 

approximation  is  valid  increases,  as  expected,  with  the  order  of  the 

network.  It  can  be  seen  from  Figures  7  and  8  that  the  approximation  in 

the  real  and  imaginary  parts  of  z ^  ,  the  input  impedance  of  the  antenna, 

is  not  as  good  as  that  obtained  for  Y  .  This  can  be  explained  by  noting 

that  the  error  in  approximating  Im[Y  ]  is  not  controlled  in  the  present 

procedure  and,  therefore,  it  contributes  to  the  observed  error  in  Z. 

in 

when  transferred  through  the  transmission  line.  To  avoid  this  effect, 
it  would  be  necessary  to  control  both  the  real  and  imaginary  parts  of 
Y  .  This  could  be  accomplished  by  using  the  ladder  element  values  obtained 
here  as  initial  values  in  a  computer-aided  design  procedure.  In  this 
event  a  nonlinear  function  minimization  algorithm  would  be  required. 

In  conclusion,  it  has  been  shown  that  the  biconical  antenna  can  be 
effectively  modeled  by  a  transmission  line  terminated  in  an  LC  ladder 
network  with  a  resistive  load.  It  has  been  shown  that  the  modeling 
problem,  can  be  reduced  to  a  straightforward  numerical  approximation  pro¬ 
cedure  followed  by  a  direct  synthesis  algorithm. 

3 .  Approximation  by  Rational  Functions 
3 . i  Theory 

It  is  often  the  case  that  one  would  like  to  express  the  transfer 
function  of  a  linear  system  as  the  ratio  of  two  polynomials.  This  form 
is  preferred  as  it  lends  itself  to  linear  transform  methods  of  solution. 

Of  the  techniques  that  have  been  developed  to  fit  experimental  data  by 
such  rational  functions  [  3 ]  —  [  4 ]  ,  the  one  by  Levy  is  the  most  notable 
1 5 ] —  [  G  i  and  forms  the  basis  of  the  rational  approximation  method  examined 
in  this  report. 
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A  fur 
H  (  j  ..  ) 


o ; 


i e  to: 


ao  + 


b0  +  +  ■■■  + 


is  chosen  to  approximate  (in  the  leas 
set  o: 


_  ■  N  (  j  ,  ) 
m  D  (  j  „■  J 


lares  sense)  a  Given  complex 


lata  F.  =  F  .  +  jF  i=l,  ...  ,  N,  where  H(j.j)  represents,  for 

1  Ri  +1 

;xamole,  the  transfer  function  of  a  lumped  network  and  F.  the  sreadv- 

i 


st; 


:ata 


scattering  object. 


,  for  example,  with  the  current  at  some  point  on  a 
coefficients  are  found  bv  mininizina 


The  a .  an 
i 

.  9 


V'. 


V 


-  H  ( j  z  i )  1  “  =  ■■ 
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.  T 


1  ' 


i=l 

The  problem,  with  this  formulation  is  two- fold:  E  is  a  nonlinear  fur.c 
of  the  unknown  coefficients  and  the  low  frequency  data  is  not  weighted 
sufficiently.  As  a  result,  wide  swings  in  the  input  data  will  cause 
large  approximating  errors  at  low  frequencies.  These  problems  may  be 
remedied  bv  defining  a  new  error, 


e*  _  D"  {  j  ,v  )  _ 

.  =  — — — — ■ - - e  .  , 

i  —  i  •} 

D  *(]J 

where  the  superscript  k  refers  to  the  iteration  number.  If,  after  each 

iteration,  one  refines  the  error  estimate  in  this  way  and  minimizes  again, 

a  much  better  approximation  is  obtained.  Sanathanan  and  Koerner  [7] 

v  v  - 1 

have  shown  that  D"  =  D'  after  a  sufficient  number  of  linear  iterations. 
With  this  change,  the  object  function  now  becomes 

'j  ;  [°R(j^i)  +  j  D  j  ( w  i  }  ]  [  F  R  (  j  w  i  )  +  j  F  x  ( w  i  )  ]  -  [  N  R  (  j  w  i  )  +  j  N  T  ( w  i  )  ]  | 

^  a  _  y  |  —  '  ~  ’  ~  ‘  —  v- 1  ‘  ~  "  I 

i-  1  D ” ‘  (  j  w ^  ) 

=  J,  I  +  jDj]  tFR  +  j  F  x  ]  -  [Nr  +  j  M  -j- 1  j2  Wik  , 

•where  W  -  •• - —  is  a  weight  function,  and  the  subscripts  R  and  I  indi- 

:  !b->i2 

male  the  real  and  imaginary  parts  of  the  terms.  The  minimization  of  E* 
at  each  iteration  is  now  a  linear  problem.  To  this  end  E*  is  partially 
differentiated  with  respect  to  each  of  the  polynomial  coefficients  and 
equated  to  zero.  This  yields  the  following  matrix  equation: 
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R,  and  I,  are  the  real  and  imaginary  parts  of  the  transfer  function  at 
experimental  points,  and  L  is  the  iteration  number.  The  coefficients 
b^,  b^»  •• •  evaluated  at  the  L-l  iteration  are  used  to  refine  the 

weighting  function  WT  for  the  next  iteration. 

A  FORTRAN  program  has  been  written  implementing  the  above  complex- 
curve  fitting  algorithm. 

3 . 2  Applications 

The  aforementioned  method  was  applied  to  the  data  generated  by 
the  terminating  admittance  function  Y  of  the  biconical-antenna  model 
as  given  by  Tai  [1].  The  data  is  shown  in  Figure  2.  A  rational  func¬ 
tion  witn  eighth  order  numerator  and  ninth-order  denominator  was  chosen 
to  fit  this  data  over  a  range  of  normalized  frequency  0  <  <  15.  This 

choice  of  transfer  function  was  based  upon  the  results  of  tests  using 
the  same  program  to  fit  the  input  impedance  of  an  ideal  transmission- 1 ine 
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termi.nated  in  a  resistance.  The  results  of  the  approximation  are 
shown  in  Figure  9  and  Figure  10. 

Another  test  was  performed  on  data  describing  the  input  impedance 

of  the  biconical  antenna  as  shown  in  Figure  11.  This  example  was 

approximated  by  a  ninth-order  numerator  and  tenth-order  denominator. 

The  results  are  presented  in  Figure  12  and  Figure  13.  It  is  seen  by 

comparing  these  results  with  Figure  11  that  the  approximated  imaginary 

part  of  Z.  fails  to  fit  the  data  near  3  2  =  0.  This  is  due  to  a  sole  at 
m 

zero  which  the  data  contained  that  the  rational  approximating  function 
could  not  accomodate  due  to  its  chosen  structure.  This  could  have  been 
corrected  by  changing  the  data  and  reinserting  the  pole  later,  a  tech¬ 
nique  described  by  Levy  [5]. 

The  poles  and  zeros  of  the  approximating  function  were  extracted 
by  standard  techniques  and  compared  with  those  found  by  Tai  and  Cho  via  a 
grid  search.  The  result  is  shown  in  Figure  14.  The  poles  and  zeros 
reflect  the  closeness  of  the  fit  over  the  approximating  range. 

3 . 3  Conclusion 

The  rational  function  approximation  method  employed  here  has  the 
advantage  of  being  able  to  produce  an  analytic  representation  of  data 
that  is  amenable  to  linear  transform  methods  of  solution.  Furthermore, 
the  implementation  of  the  method  is  straightforward  and  computationally 
efficient.  For  the  accuracy  achieved  here,  the  typical  run  took  2  CPU 
seconds  (Amdahl  470)  and  cost  $.50. 

As  a  final  remark,  a  more  nearly  mini-max  approximation  could  be 
obtained  by  incorporating  Lawson's  algorithm  [2]  in  the  iterative  pro¬ 
cedure  but  our  investigations  to  date  do  not  indicate  that  this  will  be 
necessary . 
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ABSTRACT 


This  report  deals  with  the  surface  fields  or  the  surface  charge  and 
current  densities  on  a  right-angle  earner  reflector  induced  by  a  polarized 
uniform  plane  wave.  Equivalent  circuit  parameters  of  a  short  monopole 
and  a  small  semi-loop  mounted  on  the  wedge  are  derived  and  explicit 
correlations  between  measureable  quantities  and  local  surface  fields 
are  established. 
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1.  INTRODUCTION 

The  class  of  problems  related  to  the  measurements  of  the  electromagnetic 
field  quantities  have  attracted  the  attention  of  many  engineers  engaged  in 
sensor  research.  The  major  difficulty  encountered  in  the  measuring  process 
of  any  energy  related  physical  quantity  is  the  interaction  of  the  measuring 
device(s)  with  the  physical  field  that  always  produces  a  perturbation  of 
the  field.  Therefore  it  is  essential  to  have  an  apriori  estimate  on  the 
amount  of  the  extracted  energy  by  the  sensor  and  the  extent  of  the  pertur¬ 
bation.  In  this  report  we  will  consider  a  cannonical  problem  of  this  class, 
namely  the  problem  of  measurement  of  surface  fields  on  a  right-angle  corner 
reflector.  A  good  account  of  the  previous  work  on  closely  related  problems 
is  given  in  [2], 

The  geometry  of  the  problem  is  shown  in  Figure  1.1.  The  walls  of  the 
wedge  are  perfectly  conducting  and  the  medium  of  propagation  of  the  waves 
is  air  with  oarameters  (z,  u,  a  =  0).  We  will  assume  that  the  illuminating 
polarized  uniform  plane  wave  is  propagating  in  a  plane  normal  to  the  axis 
of  the  wedge.  Furthermore  only  linearly  polarized  waves  will  be  considered 
with  polarization  of  the  E-field  perpendicular  and  parallel  to  the  axis  of 
the  wedge  resoecti  vel 


Figure  1.1.  Plane  wave  illumination  of  a  corner  reflector. 
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2.  SURFACE  CURRENT  AND  CHARGE  DENSITIES  FOR  POLARIZATION  PERPENDICULAR 
TO  THE  AXIS  OF  THE  WEDGE 

To  obtain  expressions  for  surface  current  density  K  and  surface 
charge  density  j  let  us  replace  the  above  problem  with  an  equivalent 
problem  as  shown  in  Figure  2.1.  The  time  dependence  e'^  will  be  under¬ 
stood  throughout  the  report.  Then: 

E1  (R,t)  =  Re  E1(R)ej'Jjt  =  Re  E^eJUt  -  I  •  R)  (2.!) 


X 
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l<2  =  3  ini  x  -  ccs 6  z 

=  sin$  x  +  cosi  Z 

k^  =  -sind  x  +  cos  0  z 

E  1  =  E  ^  {-cosi  x  +  s ini  z) 

o  o  ' 


(2.2) 


E2o 

■  E2o(-.OT9  X  _ 

sind  z) 

m| 

GJ 

O 

-  E3o(™9  *  “ 

sini  z) 

^  1  n 

=  c  X  + 

'40  ' 

rf  z) 

T  <J 

For  surface 

current  density  we 

have: 

K  = 

"  x  "total  *  n  x 

-(k  x  E1 
n 

where : 

n  = 

AC  3 

x  ^2  T  1^2  A  l3  4  4J 


is  the  intrinsic  impedance  of  air.  On  the  surface  x  >  0,  z  =  0  we  have: 


K(R>  =  z  x  i(k  x  E0'  +  k4  x  e4o)eJkX +  (k,  x  E,„  +  k,  x  E,„)e'jkx3i’13' 


4  2H  f(ejkX3in6  t  e-Jkx3iK9) 
0 


'2  2o  3  3o' 


>(2.3) 


■x  4H  V?ra(k xsini) 
o 


on  the  surface  x  =  0,  z  >  0  we  have: 


K(R)  =  i  x  ^  Ho1(2ejkz''"' 


39  +  2e_jkz-y0) 


z  l\'t\\jcs{Y.zcosi) 


(2.4) 


To  obtain  the  surface  charge  density  a  we  will  make  use  of  a  boundary 
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condition  derived  from  the  continuity  equation.  When  one  of  the  media 
is  a  perfectly  conducting  surface  this  boundary  condition  reads: 


v  •  K  = 

Therefore,  on  the  surface  x  >  0,  z  =  0,  we  have: 

c(R)  =  i  4H  Yjfuv  3in(kxjfna) 

*)  o 

j(R)  =  j4eE  sfK(kxjfKd) 

similarly  on  the  surface  x=  0,  z  >  0  we  have: 

-(R)  =  -j4eE  Vijz'sQ  sin(kzoos$) 


(2.5) 


(2.6) 


2.1.  Surface  current  and  charge  densities  for  polarizations  parallel 
to  the  wedqe  axis 

- — —  ,  _  -  —  — mt  — — — . 

For  this  polarization  as  it  is  apparent  from  Figure  2.2  the  following 


Figure  2.2.  Equivalent  problem  obtained  by  using  image  field. 
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H  1 

0 

=  x  -  jzni  z) 

H0 

Zo 

=  Hq"1  ( -jcj'b  X  -  sir.  )  z) 

H- 

3o 

=  H  "*  (-•_.?<> e  x  +  z) 

O 

=  H  ^  ( -jcr>~  X  +  Ikv  z) 

K(R) 

=  yj4HoV’c.;-3 

X  >  0,  2  =  0 

(2.7) 

K(R) 

=  yj4HQ  szm  szk{  kZJ.;av) 

x  =  0,  z  >  0 

(2.8) 

j(R) 

=  0  for  ei ther  x  >  0,  z  = 

0  or  x  :  0,  2  >  0 

(2.9) 

3.  OPE.’I  CIRCUIT  VOLTAGE  OF  A  SHORT  MOriOPOLE  MOUNTED  ON  THE  RIGHT- 
ANGLE  CORNER  REFLECTOR 

The  geometry  of  the  problem  and  the  significant  parameters  are  shown 
in  Figure  3.1  We  assume  that  .the  p 'obe  is  electrically  small.  We  will 


ngure  3.1.  Short  monopole  mounted  on  a  wedge. 


re¬ 


consider  only  the  non-trivial  polarization  of  the  incident  field.  As 
a  result  of  the  image  theorem,  the  induced  current  in  the  receiving  antenna 
is  related  to  the  open  circuit  voltage  by: 

V  r 

Ir  =  -r  °^-  - -  (3.1) 

Z  /2  +  Z. 
in  L 

where : 

r 

I  =  current  at  the  base  of  the  receiving  antenna 

Z.  =  input  impedance  of  two  transmitting  parallel  dipole  antennas 

resulted  by  removing  both  conducting  planes  and  using 
the  image  probes. 

r 

VQc  =  open  circuit  voltage  of  the  receiving  monopole 

In  this  section  we  will  deal  with  open  circuit  voltaae  only.  The 
inout  impedance  problem  will  be  discussed  in  a  separate  section.  Using 
the  vector  effective  height  h^  of  an  antenna,  we  can  write: 

Vq£  =  F  •  h1  (3.2) 

The  vector  effective  height  o^  a  short  dipole  (or  equivalently  a  inonopole 
on  a  ground  plane)  is  given  by: 

hfc(R)  =  .sin 0  0 


The  pertinant  polarization  of  the  incident  field  for  this  configuration 


of  probe  is  dern'cted  in  Figure  3.3.  Following  the  same  procedures  of  the 
previous  section  and  considering  the  direction  of  the  current  density  of  probe 
we  can  write: 


R‘  =  d,  9'  =  J,  y  =  0) 

\ljc  =  -  2^(y~a2)  <:  (kdo m:9)  (3.5) 

4.  INPUT  IMPEDANCE  OF  THE  PROBES 

Equation  (3.1)  of  section  3.  indicates  that  for  completion  of  the  equi¬ 
valent  circuit  parameters  we  need  to  evaluate  Z1  .  which  will  be  simply 
referred  as  Z.  from  now  on,  in  each  of  the  probing  configurations  considered 
previously.  Let  us  once  more  recall  that  Z-n  is  the  input  impedance  of  the 
transmitting  antenna  in  the  presence  of  its  images.  With  this  in  mind  we 
will  begin  deriving  analytical  formulas  for  the  input  impedance  functions 
involved  in  the  problems  at  hand  under  the  assumptions  imposed  on  electrical 
sizes  of  the  probes.  The  problem  of  determining  the  impedance  for  the  cases 
considered  here  has  been  extensively  explored  previously,  however  most  often 
in  the  form  of  tables  and  curves.  We  will  include  here  the  complete  express¬ 
ions  for  these  functions. 

4.1.  Impedance  parameters  of  two  identical,  parallel,  and  short  transmitting 
dipoles. 

The  problem  arising  from  the  application  of  image  theory  to  two  antisymmetr 
ally  driven  antenna  is  shown  in  figure  4.1.  Since  the  probe  is  assumed  to  be 
thin  and  short  (a  k l  1),  a  linear  current  distribution  is  a  suitable 


v  r  = 

oc  ^ 


~i  ht: 
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aoorox imati on  and  the  induced  EMF  method  can  be  applied  successfully  to 


determine  its  input  impedance. 

For  convenience  we  will  proceed  by 
assuming  a  sinusoidal  current  distribution. 
At  the  final  stages  the  results  wi i i  be  sim¬ 
plified  by  using  the  conditions  imposed  on 
a,  ,  and  o. 

Based  on  filamentary  current  distri¬ 
bution  we  have  for  the  magnetic  vector 
potential  A(R): 


Two  identical , parallel  transmitting 
antennas. 


R-|  and  are  the  instances  from  the  observation  point  to  the  source  points 
on  the  axes  of  the  dipoles.  For  R  on  the  surface  of  one  of  the  antennas  we 
ha  !/e : 


Rj  =  /[(z  -  z1)2  +  a2]  R^  -  /[(z  -  z’)2  +  o2]  for  p  >>  a.  (4.1) 


The  electric  field  on  the  surface  of  right  hand  side  antenna  is  given  by: 


E  -  -  jv(l  +  '~v  7-j  A  = 
,  2 


-  z  Vo 


-  Z  jrj(l  +  -2  A  (z, 

k  3z 


(4.2) 


where  as  usual  a  slice  generator  has  been  assumed.  Multiplying  (4.2)  by  I(,z) 


-li¬ 


ana  integrating  over  the  source  region  we  obtain 

vo  =  if/X1(z)I(z')!l  +  IzrfHs,  -  g2)  dzdz' 

-  a,  -  .4  K  o  Z 

-jkR-j  -jkK2 

where  G^z.z')  =  —  G.,(z,z'l  =  — 

G(z,z'  )  =  G1 (z ,z* )  -  G2(z,z' ) 

The  first  integration  with  resoect  to  z  amounts  for  evaluation  of  z-component 
of  the  near  zone  electric  field  due  to  current  distribution  of  the  form: 

I (z)  =  I  3ink[Z  -  jz: )  \z\  i  1 

m 

Therefore  it  is  given  by  [5.1]: 

l  2  I 

/  l(2l(l  +  ~~n  ~~  5" )  G  ( z ,  z  ' )  dz  =  {G(;.,z')  + 

k  K 

+  G(  -  l,z')  -  loockl  G( 0,z 1 ) } 
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Let  us  note  that  from  the  circuit  relations  for  these  antennas: 


V1  Z11  !1  +  Z 1 2 1 2 


V 2  Z21 1 1  +  Z  22 1 2 


It  follows  that  for  antisymmetrical ly  driven  antennas  we  have: 


Zin  "  Z11  "  Z1 2 


(4.4) 


Let  us  define: 


c  ^ 

!a)  =  J  {^U.z)  +  G,(-  Z,z)  -  2 -3k?.  G .  ( 0 ,  z ) } 


X  sink {Z  -  z)  dz 


(4.5) 


combining  (4.3)  -  (4.5)  we  obtain: 


-11 


Jn  r 
2 

Iv.sin  ki 


(a) 


-12 


2 

2tt3zk  kl 


Ao) 


(4.6) 


Therefore  the  problem  of  input  impedance  reduces  to  the  evaluation  of  c(a), 
and  ; ( o ) . 

t(a)  has  been  previously  evaluated  by  expressing  it  in  terms  of  sine  and 
cosine-integrals  [  6]  We  will  choose  another  approach  which  enables  us  to 
obtain  a  series  expansion  for  ?(a),  and  rjp)  and  in  particular  to  reduce  the 
results  to  simplified  forms  under  certain  assumptions  on  the  parameters.  We 


have : 


-13- 


f  C^(:.,z)  sink  (Z  -  z)dz  =  J 


-jk/(x2  +  a2) 


o  /(x2  +  a2) 


/  S-j  ( -  i,z)sink{l  -  z)dz  =  f 


81  -jk>/,(x2  +  a2) 


l  v'(x2  +  a2) 


f  r1}(0,z)sink(i  -  z )dz  =  fQ 


-jk/{x2  +  a2) 


0  /(x2  +  a2) 


Therefore  we  are  led  to  define: 


l  e-k/(x2  +  a2) 

S  (  ,  V*  )  ~  J  ry  ry  S  U'i  k  X  d  X 


*o  /(x2  +  a2) 


f 


0  e-jkAx2  +  a2) 


o  /(x2  +  a2) 


;kxdx 


Then  in  terms  of  sU,a)  and  C(2,a)  we  can  write: 


sin kxdx 


sink {2?.  -  x)dx 


r. k(:c  -  x)dx 


(4.7) 

(4.8) 


;(a)  =  s ( Z , a )  +  si^2k2(C(2T,a)  -  C ( Z ,a ) )  + 

-  ?ss2kZ(s{2 1, a)  -  s(?,,a))  -  2sinkZ3oskZC{Z,a)  +  2  cos^kl  s ( 7 , a ) 

■; ( a )  =  2  sU.a)  +  (2  s(t,a)  -  s(2.?.,a ))sos  2k Z  -  (2CU,a)  -  C{2l,a)s  ;-:2kZ 

(4.9) 

In  order  to  evaluate  s(?-,a)  and  c(t,a)  let  us  introduce  the  following  dimension¬ 
less  parameters: 


%  =  ki 

.  Gt  =  j  ,  0  =  Al  +  r 

x2) 

(4.10) 

Then : 

s  (■•  ,a  ) 

J  e-.i-;/(x2  +  a2) 

sinfxdx  =  s(  -;vv) 

(4. 7') 

o  /(x2  t  i) 

c(r-  ,a  ) 

_  f  e-M*2  *  *2> 

en.nfxdx  =  c  (f,a) 

(4.8') 

o  /(x2  +  a2) 
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wi  th  : 


3s 

a ' 


J  e-j-'^2  +  «*). 
0 


"xax 


Jxe-i-yf^2  +  A  . 

J  — • — — - xax 


'o  ^x2  +  a2) 


1  xe-iC/(x2  +  /) 

J  - 5 - -2- - ^nxdx 

o  lx"  +  a) 


1  -  j  ^  x 2  +  a2) 

—  e  J  ’  1  cos-rx 


1  e-j-y(x2  +  a2)  sin^xdx 

J  o 


Hence: 

||=  ^  +  **>  cos' x 

os  J-, 


J  “3 


sirni 1  a  r 1 y : 


■3C  r 

H  "  ‘  J  /  0 


1  -^/(x2  +  "2)  txdx  -  / 


1  xe-J7<x2  +  a2) 

- 2 - 2~ —  x  sin^xdx 


o  /  (x  +  ct  ) 


ac 


H  j-; 


e-i--  ^  +  At*. 


I_  p-j5S,xr; 

J 


(4.12) 


p  ^  /*• 

We  will  proceed  by  finding  a  series  expansion  for  —  and  For  the  sake  of 

numerical  comoution  we  will  develop  two  different  series  expansion  depending 
on  the  relative  values  of  at  =  ka. 


|1  =  _  1__  ;  L  (e-js(3-l)  +  e-j-;(3+l)j  _  e-j:,  t, 

os  Js  2 


Z  {(3  -  l)n  +  (3  +  l)n  -  2,n} 

2  n=l  n! 


s(:,v) 


W  ?  <  (S  -  U"  +  (3  +  l)n  -  2an: 

J  n=l 


(4.13) 


-15- 


.  s  (  2 r  “)  -  -  J-  ?  LdUll  s 

S  ^ ’  2  j  2  rin  n 


•n  =  2 1  '  n { ( 3- 1  )  n  +  ( 3+ 1  )  0  }  -  {(3(|)-l)n  +  (B(|)  +  l)n} 


0 1  -  n  n 
-2  a 


(4.14) 


Similarly  we  have: 


±c  =  p- 3 

3  C  j  v 


e'j"3^nC  =  J  I  1  U  S  - 1  )  °  -  (3  +  1  )  °  } 

n  =  1 


,  ®  i  .r,n 

c(;,-i)  -  c(o.a)  =  -i  i  -jfjr  -  ( e+i )  > 


c  (  0  , 


/  - ^—9—  =  c  (x  +  /(x2  +  a2)) 

5  /( x2+a2 )  n  0 


.  co  /  .  \  rl 

etc,,)  .  ;.n  iii  -  1  l  ((3-1)"  -  (p-M)"} 

a  n=i 


(4.15) 


a  2[3(f)  +  l]  -i  «  (-i2~)n 

2  c  (  '  ,  :t )  -  c  (  2  -  ,  f)  =  J,n - o~  +  £nct  -  ■=■  £  — 


(3  +  D 


n  1  n  n 


(4.16) 


c  =  21'n[(3-l)n  -  ( B+ 1  ) n  ]  -  [(3(J)-l)n  -  (3(|)  +  Dnl 


The  previous  expressions  are  so  far  exact.  Let  us  introduce 
aDproxi mat  ions  under  the  assumptions: 


1  3  =  kZ  <<  1 


then 


3(a)  =  1  +  2""  1 2  +  0  ( ci  ^ ) 


2(2(5-)  +  !] 


(  3  ( a )  + 1  ) 


- o  =  in  4(1  +  fa  +  ...)  -  2 £ n 2 ( 1  +  t-  +  ...)  = 


0  (  a  2  ) 


I  (")  6m[l+(-l  -  I  (")  b'"(£)x 


n)  Bm(f' 

m  2 


n  x  t  n  o1  'n  n 

x [1  +  ( -1 )  ] -2  a 
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„  1  -  n  m  , 

L  3  l  Jt 


-  Sm(|)  =  2 1  ~ n  (  1  +  a2)  2  -  (1  +  \)2 

=  ( 21 "n  -  1 )  +  0(a2) 


3  =  "  (21_n  -  1  )  (  n )  [  1  +  (-l)n"m]  -  ctS  ,  +  0(ct2) 

n  m  nl  v  ' 

m  =  0 

=  ( 2 1  " n  -  1)  2n  +  0 ( a2 )  =  (2  -  2n)-a6nl  +  0(a2) 
where  5  ■]  is  the  Kronecker  delta. 


c  =  -2 
n 


l  (")  3m[l  -  (-I)"""1]  +  l  (")  Bni(|}[l  -  (- 


n  \  0  m  /  ol  • 
m>.  3  (2 


1  )n_rn] 


=  (2n  -  2)  +  0 ( a 2 )  . 

Thus  we  have  the  following  approximations: 


*(:.=•>  •  4  (  I  +  +  0(a2O 


(  4  .  1  3  1  ) 


2  s  (  ,  a )  -  s  (  2  T  ,  f)  =  {  I 


&  \  _  1  r  v  (  ~  1  2  c  )  /0n 


(2"  -  2)  -  j  2 a } 


(4. 14’  ) 


+  0(a  >) 


(4.14’  ) 


2  c  (  ,  a )  -  c(2^,  f)  =  2na  -  1  £  (~^2^-  -  (2n  -  2) 


+  O(oc^)  . 


Finally  for  t ( a  )  we  obtain: 


°°  l  •  p  ■)  n 

'(a)  =  3i>:Z%Zna  +  ( sinZZ  ~  joos2r)  £  -  x 

n  =  l 

x  (2n  ^-1)  +j  l  - — ^  2-r-^—  -  (2  +  cos2')aC  +  0(a2- 
n=  1 

t(a)  =  2 r  -  ^'(2  0 3  -  3a;  +  3i>i2Z,lna  -  ^(2l)4  + 


( 4. 1 6  1  ) 


+  0  [nax{  -J*  ,a2>  )  } 


(4.17) 


>>  1 
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For  this  case  we  have  %  =  k  l ,  ct  =  - .  Following  the  above  formulations 


we  set: 


s  (  v  ,  :i )  =  s  (  0  ,  at )  +  e 


-J 


I  8»5 

n  =  0 


(4.18) 


C  (  :,  5  ,i  )  C  \  0  ,  Cl  )  +  e 


-  J  \a 


i  c 
n  =  0  ' 


(4.19) 


Then: 


is  - /)  -  j.  v  _ n -  1 

i  -  L ,  n 

’  n  =  1 


-  jae"J-a  I  s  rn 
n  =  0  n 


-  J 


i  [-j;l3  <-  (n  +  l  )«n  +  1  ) 

n  =  0 


_  n 


=  _  _ L  (e-^3- 

j  v  e 


.  e-j^  -  ^  x 

6  '  323  nil  nl 


x  { ( 3  -  a  -  1  ) n  +  ( 3  -  a  +  1  ) n } 


n 


°°  ,  .  x  n  - 1 

=  p  ~  -1  ’ "t  V  (  ~  n  =  p  J  i  -i — -  n 

n  =  1  n  ■  n  n  =  0  rn^tlT  n  + 1 


r  (  -  J  O 

L  T~  v'-'  a 


where  we  have  defined: 
1 


'  n  2  v  “ 
Therefore  we  have: 


ct 


-  1  )"  +  (3  -  a  +  1  )"};  n=  0 , 1  , 2 , 


(4.20) 


■  jots  +  ( n  +  l )s  ,  , 
J  n  v  '  n+l 


=  (-j)n  ^nTTTTi  n  =  0 , 1  , 2 , 3  , .  .  . 


(4.21  ) 


Recurrsion  relation  (4.21)  starts  off  from  Sq  =  0  which  is  an  immediate 
consequence  of  (4.18): 

-n  =  0 


1  =  "l 


it  2 


2!  '  2  1 


'J2 
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s  =  HU  (Zl  .  ;t_^2  +  2  ) 

3  3 !  1  3  2  ~  r 


■  j  )  n ~  (£n  _  t?n-1  :t  ~n  -  2 

n  -  2 


n  ! 


n  n 


,  r  i  \  n- 1  n  - 1  » 

+  (  -  U  a  a1  ) 


n  =  1  ,2,3,.  .  . 

3  c 

Similarly  for  —  we  have: 

d  ' 


(4.22) 


-  l  -  a 


n  =  0 


+  (n+l)c„+1) 


=  L.  e'j 

J  s 


i  •  -  00  t  ■  -  \ n 

r  =  _  I_  p-J-Ci  c  (  -  J  -J  r 


2-: 


n  =  1 


n  I 


(2  -  a  -  1) 


-  (e  -  a  +  1 )n] 


00  /  -  -  \  n 

=  e  J  £  -r — -V-v-:-  Y 


n  =  0 


( n  +  1  )  !  f  n  + 1 


(4.23) 


w  i  t  h 


V  n  "  { ( 3  ~  ct  +  1  )  -  (3  -  a  -  1  ) 


(4.23) 


Finally: 


cn  +  (n  +  1)  cn  +  1  =  ( n|j  )  ,  Yn  +  1;  n-0,1,2... 


(4.24) 


C0  =  0 


C1  =  Y1 


1  Y  ?  aY 

.  (_ £  _  - L) 

2 1  v  2  i  ' 


1 


/  .  N  n  -  1  Y  a  Y  , 

_  (  -  J  )  t  n  n  -  1 


n  -  1  n  -  1 


n  J 


<-£  -  -sir1  +  ■■■•  +  <-’>  “  '' 


(4.25) 


In  general  (4.18)  -  (4.20)  together  with  (4.22),  (4.23),  and 


(4.25)  are  compact  enough  to  permit  numerical  computation  of  s ( £  ,a ) 
and  c(5,t)-  However  we  would  like  to  discuss  the  approximation  and 
further  simplification  under  the  conditions: 
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:t  =  •  <<  1  ,  ■  <<  1  • 

We  have: 

2,  _  v  /  ,  x  m+ 1  (  2m  -  3 )  !  !  .-(2m-l) 

b-a  =  v'(l+a)-a  =  )  (-1)  - - - -  j. 

m  =  1  2mm I 


with: 

(  2  :n  -  3  )  !  I  =  (  2m- 3  )  (  2m- 5  )  -  -  -  3.1  m=2,3,... 

( -1 ) ! !  =  1 
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1 

J 


Y  2 

_  2,  . 

■  K 

1 

4 

-3  .  1  -5  5  -  7 .  ( 

Y-j  ’ 
J 

-r[l  +  3 
J 

{  = -2t)  ~  ) 

4'1 

.  3  -2  3-4.  ,  -6. 

+  4  x  -  g  X  )  +  0  ( u  } 

v  a 

x)  +  4(2 

'  3 1 

-  Li  J  ] 

=  j(2-t'1)  +  OU'5) 

Fro 

n  (4.22) 

,  (4.25) 

and 

(4.26)  it  follows  that: 

/ 

'  ’  1  ' 

. )  =  i  2  x 

r1  -  (2 

-“i 

- .  N  ”  4 

1 

+  0(,'5) 

’  ?  v 

\  -  j 

'  '  2 

7  [(2:0' 

2  -  2  (  2  a  )  '  4  ]  +  0  (  a  ~  5  ) 

■t( 

J 

\  _  1 
■r  '  '  '  3 

7  [  2“  \  2  ,  j 

-1  + 

3  (  2  .t )  “  3  }  +  0  ( :c "  5  ) 

-•4:  U  =  j}  [f(2:0"2  -  2(2, )'4]  +  Of,'5)  (4.27) 


C  „  ( 
c 


c ..  (  ( ;  = 


c , 


=  -  ~(~x  +  -  ( 2  ~ 3 )  +  0(a'5) 

-  fp'  '.2  -  ~  +  f(2  x)'2  -  3(2, )"4  +  0(a"5) 
|t ^  3  +  f  +  5-(2:0^  +  f  (  2  :t )  "  3  +  0(:t"5 


.  1 s  c  ,  let  u s  note  that: 

1  CO 

c  (  0  ,  i  +  /  ( 1  +  -X  2  ) )  -  "  n  (  1  +  j i  +  V  (  -  1  ) lT1  +  ^  x 

m=  1 


■  2 -  3  (  !  !  -  2m  x 

x  - - -  -  t  ) 

?.  ”i ! 


+  /. 
m=  1 


!-})n-]  -t-2n  _ 
'  •  m 


2  mi  ! 


r 
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N  o  w 


1,-1,  r  /  n  ,n  +  2  (  2  in  -  3  )  !  !  -2m,  2 

-  2  ■'  f  2  f-’)  ~ ~  11  '  + 

m  =  1  2m! 

-1  1  -3  .  n,  -5, 

-  -  -6J  +  0'J  > 

2c(0,:,)  -  c  ( 0 ,  |)  =  j'3  +  0!„,"s) 

we  can  proceed  to  evaluate 

=  2  s  (  •; ,  a )  +  .-s2-;(2s(;,a)  -  s(2^,  |)  ) 

-  2-  (2c  -  c(2,r,  |))  = 

W  -1 

-  J ,  i  '  C 


=  e 


n  =  0 

x  (  2  1  n  c  ( 


(  2  " )  n  {  2  1  “  n  <3  n  (  :c )  +  -jos  2  E,  x 


„•>>  ‘  VO  -  ^£«ZS(21_nc1 


(4.28) 


-  c  {%)))';  -  cin  2-;(2c(0,a)  -  c(0,£)) 


=  e'j  ^  l  l  (2-;)  '-2 

n  =  0 


n  -„2-  n 


~  3n{i]  + 


(20  2  (01)  4 

- 0“  +  24  •  - 


-)(2,-",n(c,)  -  ,n(f)) 


-  (2- 


(20 


-  O  (  2 1  ~ n  cn(-0  -  cn(|})}} 


•>  2  '  (2c  (0,0  -  c  (  0 , 0  ) 


-  J  r. 


1  4 


-  3  (  o  -  \  •  3  -2/^_\2  /  9  -  \  3  (  1  j_ 

2  (  2  : J  +  j  g  a  (  2-,)  +  (  2  :J  (  Yg  A  ;  + 


\  \ 
J  j 


(4.29) 


-  0  3  (  2  )  f  0(-;,:(  O3-  3,,050) 


f  r. 


2  i' 

o- 


,1  r1  »  ,i  f  ’  *■  f 
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'(■)  =  ^2  e  ~ j  k  (  2  £  )  4  {  (  k'  )  ”  1  +  j3(k.)"2  +  3(k.)'3} 

1  '  o  -  \  4  /  ,  \  -  3  «!-  -5  -3  ”3i  i 

-  g\2u)  (k  )  +  a  ); 


(4.30) 


using  (4.6)  together  with  (4.17)  and  (4.30)  we  finally  obtain  expres¬ 
sions  for  z  and  z ^  ^  : 


1  i 


2  r4 

2  _  ’  ( a  )  ~  2^~  ’  ^  5  +  24~  + 


)  ’ 2  x  :  (  a 


(  r2  +  j  +  OU4))  '(a) 


'1 1 


^  {(4  -  6x)(27)-1  -  (|  +  f)(2T)  + 


2  a  s  t  '  c  n  x 


12 


9  2 

2-:)c}  + 


2,  ,  «r  ./  -3  ..2.-1  v  , 


w  i  t  h 


(4.31) 


12 


J 


2-r  in 

1  -  i  k 


",(  :)  *  jr  (T2  +  3  +  OU4));'..-) 


(  2  - )  2  { (  k  .  )  "  1  +  j  3  (  k  (; )  ’  2  +  3(k.-)“3 


l(2-:)2(k  ,)'3  +  0  {  '  4-:  !  (  k  )  "  3  ,  (k  :)'5!  : 


Then  the  input  impedance  for  the  problem  at  hand  is  given  by: 


Z  ■  -  z ,  ,  -  zno 

in  11  12 


(4.4) 

Let  us  note  that  in  (4.32),  as  in  previous  sections,  we  have 


4  .2  I_n  n  u  t  i  m  n  e  d  a  n  c  e_  of  two  coaxial  small  circular  loops  in  transmitting 

mo  d_o  . 

For  an  electrically  small  loop  we  assume  that  its  current  distri¬ 
bution  is  uniform  when  driven  by  a  localized  voltage.  The  conventional 
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in  duced  EMF  method  will  be  applied  to  determine  its  impedance.  How¬ 
ever,  the  formation  is  much  more  complicated  than  the  case  for  linear 
antennas.  Therefore  assumptions  imposed  on  the  geometry  of  the  probing 
antennas  justify  adoption  of  simpler  approximate  methods  in  ad  hoc 
bases  for  evaluation  of  self  and  mutual  impedances. 

As  can  be  seen  from  the  application  of  the  EMF  method  to  the  anten¬ 
nas  of  previous  the  previous  section  the  self  impedance  can  be  obtained 


Figure  4.2  Coaxial  transmitting  loop  antennas 


by  removing  the  image  of  one  antenna  and  evaluating  the  input  impe¬ 
dance  of  the  isolated  loop. 

The  input  resistance  of  a  constant  current  loop  may  be  simply 
evaluated  by  an  application  of  Poynting's  theorem  [9.2].  The  result 
i  s  : 


R  . 
i  n 


6  - 


(  t  (  k  a  ) 2  )  2  . 


(4.33) 
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A  rather  simple  way  of  obtaining  a  compact  formula  for  the  input 
reactance  of  the  loop  is  to  make  use  of  the  reactance  of  the  loop  based 
on  the  circuit  theory.  According  to  Reference  [9.1]  the  reactance  of 


the  loop  is  given  by: 

x  j  n  -  “(Li  +  L0>  :  “!l'i  +  aU 

where  L .  is  the  internal  inductance 

selected  mutual  inductance,  and  b  i 

loop  antenna.  This  result  compares 

of  the  formula  obtained  for  x.  by 

l  n 

discussed  in  [1].  Combining  (4.33) 
impedance  of  the  loop: 


^  -  2)}  (4.34) 

of  the  wire,  L  ^  is  the  so  called 
the  radius  of  the  wire  of  the 
very  well  with  the  leading  term 
pplication  of  the  wave  theory  as 
and  (4.34)  we  obtain  the  self 


R.  +  j  x.  -  ~(  m  (  ka  )  2  )  2  + 
in  J  i  n  6 


+  j. 


+  a[:n(34)  -  2]} 


4 . 3  Mutual  impedance  of  the  two  coaxial  loop  antennas 

As  for  the  case  of  monopole  on  a  corner  reflector,  we  will  pre¬ 
sent  a  simple  and  compact  formulation  for  °f  two  small  coaxial 

loon  antennas.  He  will  use  the  EMF  method  again  and  therefore  the 
Fresnel  field  of  a  constant  current  loop  is  needed. 

Trie  Fresnel  field  of  a  constant  current  loop  antenna  has  been  pre¬ 
viously  obtained  in  the  form  of  a  rapidly  converging  power  series  in 
[10.1],  [3],  [3],  and  [4].  We  will  obtain  another  series  expansion 

which  closely  follows  the  ones  g^ven  in  [8]  and  will  prove  more  suita- 
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i  R  -  R '  !  =  /(  R2  +  a2  -  2a  R  jsay) 
where:  jooy  =  cini  jo  3  ■'  • 

We  conclude  that  ■  is  the  angle  between  two  position  vectors  v  and  / 
with  angular  coordinates: 


2  ’ 


=  0 


2’ 


Fro:"  (  4  .  3 S  )  we  have 

'V2  +  a2  + 


((R2  +  a2)2  -  4(aR aind)2  )  Y2 


=  (R2  +  a2)T^O  ±  /(l  -  <2))}h 


2a  R.-  •••:• 

R  24-q  2 


<  1  . 


(4.39) 


Now  we  can  use  the  addition  theorem  for  the  spherical  Hankel  function 

h{Q2)(k]R  -  R1  |  )  -  h^2)(k|~  -  “ ■  |  ) 

and  the  addition  theorem  for  Legendre  polynomials: 

n 


n 


.  PnU.,c9)Pn<^0')  *  2  I  Vlieoct)  x 


m=  1 


where  - 


x  pT  oo3b  1  )aosm{  0- ) 


+  sznBainB  <  aos( y-  $  •  )  =  ao3$ 


analogous  steps  which  led  to  (4.36)  we  obtain: 


AT 


•  j  u  k  a  I  n 


)  l 


4n  +  3 


n  =  0 


TTIT  (  2n  +  l  )‘L  2n  +  l 


[PL.n(O)]2  x 


x  j2n+]  ^  kv  ’  ^  h  2  n  + 1  ^  kv^ 


F  o  1  1  o  w  i  n  a 


(4.  40'  ) 


Furthermore  we  have  [ 7  j  : 
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P  4  (  0  ) 


.-1/2 


•  '■>)] 


.  1.  ,  1 

r(2  2;  2U 

r  ( l  +  -  \\x 


Therefore: 


P2n+1<°) 


Finally  we  have: 


t \ n  + 1  ( 2n  +  l  )i 
]  22n(n.)2  ‘ 


jiika  I 


0\  r  ( 4n  +  3 ) { 2n  +  l )  r  ( 2n  )  I  ,2 


•)  I 

n  =  0 


2^n(n!  ) 


x  j2n  +  l ^  kV ^h2n  +  l ^  k 


Thus  electric  field  intensity  can  be  obtained  as: 


E(R)  -  -j..(l  + 


)  A  -  -  j  A  - 


( ka )  2fil, 


0  I 


- A Jl+3 _ r  p  1  (Oil2  X 

n+1 ) (2n+l ) 1  2n+l 1 J n 


J  2  n  + 1  (k  ’’  )h2n  +  l^kV)  ''  < 


(4.41  ) 


flote  that  in  obtaining  the  expression  (4.41)  we  have  used  the  follow¬ 
ing  relation: 


7  •  A  =  -d-4— r  It  Ac,  =  0. 

Now  we  are  in  a  position  to  obtain  the  mutual  impedance  of  the 
antennas  shown  in  Figure  4.2.  Application  of  the  reciprocity  theorem 
to  one  of  the  two  antennas,  say  antenna  2,  yeilds  [5.2]. 

Vr  =  -  -1  /  E1 • dv  ' 

I  1  V' 


where 


a  open  circuit  voltage  of  antenna  2  in  receiving  mode 
a  current  of  antenna  2  in  transmitting  mode 


a  incident  electric  field  when  antenna  2  is  removed 
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J  1  =  current  density  of  antenna  2  in  transmitting  mode. 


Noting  that: 


■21  "  I 


2  OC 


1  1  I 


2  =  0 


we  can  equivalently  write: 


•21 


rV  /  E1  .  JJR'  )  d  v 

M  2  V ' 


(4.42) 


where  : 


J  2  (  R  '  ) 


I2S(: 


,  \  S  (  R  '  -  a  ) 


R 


Using  (4.41)  in  (4.42)  we  have: 


■21 


-{.< 


(ka) 


4n  +  3 


4a 


£  (n  +  T)  (  2n+l  )  [P2n+l  1  x 
n  =  0  v 


x  J 


2  n  + 1 


(kv')i>^i,(kv)e(»*3')  5(R^:al 


x  R  1  2 .3  f  9  ’  d  R  1  d  0  '  d  ' 


■21 


(  k  a  ) 


,  2 


4n  +  3 


4a  “  (  n  + 1  )  (  2  n  + 1  T  [P2n  +  l(0)]  x 


x  J2n  +  1^;  +  l  ‘ 


2  n  + 1 


') 


( p  1  -  a  )  P  , 


V 


R  1 


:9  1  dR  '  d9  '  d  - 


w  e  nave: 


21 


(ka)"  l 


4n"J  [pL.t(O)]2  x 


2'rw/  n^Q  (  n  +  1  )  (  2n+1  )  1  2  n  + 1 


j2n  +  i  (  kv  '  ^h2n  +  l  ^  k'J  ^ 


(4.43) 


For  |  k a  !  <  <  1  and  a  <<  we  have: 


2a' 


2  a  R.:  -V  : 

‘  2  2  2  2 

R  +  a  :  +  2  a 


■<  1 


(4.  39  1  ) 
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(R2  +  a2)1'2  {hi  ±  /(I  -  v2 )  j  } 


=  (R2  +  a2}’2  {^[  1  ±  (1  -  l,2  +  0(  . 4 ) )  ]  P 


2 

=  (R2  +  a2)'2  |[1  +  0(v2)]  =  — ^ — y 

( ■'  +  2  a 


r  H  +  0(  :■  2 )  1  (4.44) 


(R2  +  a2 )  2  [1  -  k2  +  0(v4)]  =  (:2  +  2a2)’2  [1  +  0(.2)]. 


Therefore ; 

k'.  '  =  ka 


n  +  0(.2)]  -  ka  — ^ - — -  n  +  0(k)] 

.?  2  +  2  a  2  )  2  (c2  +  2a2)  2  4 


k  .  -  k  ( .'  2  +  2a2)'2  [1  +  0  (  ; 2 )  ]  =  k(,2  +  2a2)’2  [1  +  0(^j)]  (4.-4'  ) 


Let  us  note  that  [11] 


i  r7\  -  (D-vS  1  (y)  -  9n7n  V  (  -1  )m(  n+m)  1  2m 

J  n  v  z  '  ^  2  z  J  n  +  v  z  '  1-2  4  m  (2n+2m+l  )  !  Z 

n  =  0 

(  2 ) 

The  asymptotic  expansion  for  h  ^  1 { z )  for  large  argument  z  reads: 

>'2)U)  -  H - 1  e-Jt«-<n*Di)  , 


J- 


F-l  (  - n  )  (  n  + 1  ) 
l  ■  ^ . ~m~~  +  0(2 


m  =  0  m ! ( -2 j  z ) 

We  conclude  from  the  above  that  under  the  conditions  imposed  on  a  and 
the  first  term  of  the  series  (4.43)  will  be  sufficient  for  computation 
of  z, 


■21 

z 


21  =  t|i:(ka)2(  pJ  (0)  )2j1  (leu  1  )h1(2)(kv  )  + 


+  0  { (  k  a  )  2  (  kv  ‘  )  3  (  kv  )  ”  ^  } 


z21  =  n ka  )  2  j  ]  (  k;  1  )  h  \  2  }  (  kv  )  +  } 


(4.45) 
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'21 


T  /  ,  \2  >. > 1 

=  -n  o(  ka  )  e 

v 


'J'kv  +  o  { (  k  a~-/r> 

(ko)4 


f(ka)2  2 


e  x  p  { -  j  k  (  p2  +  2a6)2} 


o  +  a 


t  0(U“LL> 

(H) 


(4.46) 


Let  us  once  more  recall  that  the  impedance  which  should  be  used 
in  (3.1)  can  be  obtained  from  (4.35)  and  (4.46)  by  means  of: 


z  . 
l  n 


"11  1 2 
and  changing  c  to  2d 


1  1 


21 


4 . 4  Input  impedance  of  two  antisymmetrically  driven  identical  coplanar 

circular  loop  antennas. 

The  self  impedance  for  this  configuration  of  loops  under  assumptions 
imposed  on  the  geometry  of  the  antennas  in  section  3.1  is  identical  with 
the  input  impedance  of  a  single  loop  antenna  and  can  be  obtained  from 
(4.35).  The  mutual  impedance  can  be  obtained  in  exactly  the  same  man¬ 
ner  as  followed  in  section  4.3,  however  the  integrals  appearing  in  this 
case  are  a  little  involved  and  for  the  purpose  of  the  problem  at  hand 
it  suffices  to  use  an  asymptotic  formula  for  z  which  is  developed  in 
[10.2]  based  on  effective  heights  of  the  transmitting  antennas: 


Figure  4.4  Coplanar  loop  antennas 
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,-jkP 


2  \  P  ^1^2 


(  h  1  •  h  2  )  =  j  n  -2-x~p--  x 


x  [Jk--(ka)' 


x  (f  (ka )  3i«9,] 


:  2  ‘  lj  2 


z  -j  2  _  -j  j  n(ka)' 


-jkP 


(4.47) 


5 .  Correlation  of  the  unperturbed  surface  fields  to  the  equivalent 
circuit  parameters  of  the  probes . 

In  previous  sections  an  attempt  was  made  to  completely  describe 
the  equivalent  circuit  parameters  of  the  different  probes  mounted  on 
a  corner  ref  1 ector.  Si  nee  any  physical  measurement  performed  by  the 
probes  can  be  described  completely  in  terms  of  the  open  circuit  voltage 
and  the  input  impedance  of  the  sensors,  we  will  attempt  to  relate  the 
surface  field  quantities  in  the  absence  of  the  sensors  to  the  equivalent 
circuit  parameters  of  the  sensors,  or  in  other  words  to  the  measur¬ 
able  quantities. 

In  section  2  it  was  shown  that  for  a  plane  wave  illumination  of  the 
wedge  with  electric  field  polarized  perpendicular  to  the  wedge  axis 
surface  current  and  charge  densities  are  respectively  given  by: 


-x  4  H , 


(  kx  a  >.e  ) 


0  z  =  0 


(2-3) 


o(R)  =  j  4-  E  g  o::reof>:(kxnfn0) 


>  0  z=  0 


(2.5) 


On  the  other  hand  the  open  circuit  voltage  of  a  short  monopole  mounted 
on  the  wedge  was  found  to  be: 


j2;.E 


(3.3) 


Comparing  (2.5)  with  (3.3)  we  establish  that: 


Z  >  oc 


(5.1) 
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The  unperturbed  surface  charge  density  is  therefore  related  to 
the  open  circuit  voltage  of  the  probe  by  an  equivalent  capacitance  per 
unit  area: 


A  2  ; 


'eq 


Farad/m1 


(5.2) 


Equation  (5.1)  is  the  manifestation  of  the  electric  coupling  of  the 
monopole  probes  and  further  justifies  the  name  of  'charge  probes'  given 
to  this  kind  of  sensors. 

Similarly  the  V f  for  a  semi  loop  whose  axis  is  parallel  to  the 


wedae  axis  was  found  to  be: 


ur 

'oc 


■  j  2  4’--a2  )Hj  sea  (kxair.e  ) 


Comparing  (3.4)  with  (2.3)  we  have: 

.  a-a2 
=  ■] 


r 


oc 


(3.4) 


(5.3) 


r 

That  is  the  equivalent  inductance  relating  V  to  unperturbed  current 


density  is  given  by : 


L0q  "  j  u  (  tt  a  ^ )  Henry  m.  (5.4) 

It  is  obvious  that  the  coupling  of  the  probe  to  th.e  electromagnetic 
field  in  this  case  is  of  magnetic  type.  Let  us  finally  note  that  for 
the  plane  wave  whose  electric  field  is  polarized  parallel  to  the  wedge 
axis  the  surface  currents  and  the  corresponding  surface  magnetic  field 
can  be  detected  by  a  semiloop  sensor  whose  plane  is  parallel  to  the 
current  lines.  From  (2.7)  and  (3.5)  we  obtain  for  this  case: 


V 


r 

o  c 


(5.5) 


That  is  to  say  the  equivalent  inductance  for  this  case  is  also  given 
by  (5.4). 
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Let  us  conclude  from  the  above  results  that  as  long  as  our  sen¬ 
sors  are  electrically  small  low  frequency  elements  relating  the  open 
circuit  voltage  of  the  probes  to  the  surface  fields  only  depend  on 
the  geometrical  characteristics  of  the  probes  and  are  independent  of 
the  characteristics  of  the  source. 
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